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Abstract 

First-principles calculations have been performed for different collinear magnetic orderings in 
orthorhombic manganites, such as HoMn03 , TbMnOa and YMnC>3 , showing large GdFeC>3- 
like distortions. Our results suggest that the AFM-E type ordering, experimentally observed in 
HoMnC>3 and recently proposed from model hamiltonian studies as a potentially novel phase, is 
indeed the magnetic ground state. Its stability is strongly connected with octahedral distortions 
and points to the relevance of structural more than chemical effects. The calculated exchange 
constants, extracted from a Heisenberg model used to fit the first-principles total energies, show 
that the ferromagnetic in-plane nearest-neighbour coupling is reduced compared to less-distorted 
manganites, such as LaMn03 . In parallel, the antiferromagnetic next-nearest-neighbour coupling 
along planar Mn-O-O-Mn paths in highly-distorted manganites plays a relevant role in the stabi- 
lization of the AFM-E spin configuration. In agreement with experiments, the density of states 
shows that this phase is insulating with an indirect band-gap of ~ 0.5 eV. 

PACS numbers: 75.47.Lx,75.30.Et,71.70.-d 
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I. INTRODUCTION 



Rare-earth manganites show a fascinating variety of physical phenomena--^ ranging 
from their unusual magneto-transport properties ( leading to the well known colossal 
magnetoresistance^) to peculiar charge, spin and orbital orderings to recently discovered 
ferroelectricitj*^. It is widely accepted that the intriguing manganites physics is intimately 
related to the profound interplay between the crystal lattice and the spin, charge and or- 
bital degrees of freedom. In RMnC>3 with light rare-earth (R) cations (R= La, Pr, Nd) the 
collinear antiferromagnetic A-type (AFM-A) ordering (i.e. ferromagnetically (FM) coupled 
ac layers antiferromagnetically aligned along the b direction in the Prima setting) is favored, 
whereas manganites toward the end of the series (R= Tb, Dy, Y) exhibit a sinusoidal mag- 
netic structure defined by the propagation vector (k x ,0,0). Recently, a new AFM phase has 
been observed in late-R HoMnC>3 5 as well as in nickelates 6 , with an "up- up-down-down" 
spin ordering in the MnC>2 planes (the so called "E-type" in the Wollan-Koehler notation^), 
i.e. two Mn with up spins are alternated by two Mn with down spins along the principal 
axes of the in-plane cubic unit-cell, whereas the spins are reversed from plane to plane in 
the perpendicular direction. It is surprising that only very recently this phase has been 
proposed in undoped manganites^' 9 , considered well-understood systems, where the E-type 
ordering was found to be stabilized in a wide region of the phase space. 

The focus of the present work is therefore to compare our results with experimental 
observations suggesting AFM-E as the magnetic ground-state in HoMnOs and, more in 
general, to deeply investigate this non-conventional magnetic ordering, focusing in particular 
on the link between magnetic and structural properties. In closer detail, we consider different 
collinear magnetic alignments for three highly distorted manganites, HoMnC>3 , TbMnC>3 and 
YMnOa , in order to find the stable magnetic configuration and compare our results for the 
most studied LaMn03 . Surprisingly, we find that the AFM-E is the ground-state collinear 
phase for HoMnOs , TbMnC>3 and YMnC>3 and that its stabilization is mainly related to 
octahedral distortions that decrease the in-plane FM interaction between e g states. 
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II. STRUCTURAL AND COMPUTATIONAL DETAILS 



The simulations have been performed within the generalized gradient approximation 
(GGA) to the exchange-correlation potential in the density-functional framework. The 
all-electron full-potential linearized augmented plane-wave (FLAPW) 10 formalism in the 
FLEUR— implementation has been used. Muffin-tin radii have been set to 2.5, 2.5, 2.5, 2.5, 
2.0 and 1.5 a.u. for Tb, Y, La, Ho, Mn and O atoms, respectively, whereas the wave function 
cut-off was chosen as 3.8 a.u. -1 . Convergency was carried out using 24 special k-points in 
the orthorhombic Brillouin zone, whereas finer quantities (such as differences in the total 
energy) were checked using 192 special k-points. 

In order to accurately treat within DFT the 4/ electrons - that typically occupy an 
unphysical position slightly above the Fermi level (Ep) when considered as valence states 
- the Ho and Tb 4/ electrons were kept in the core as an open shells (with 7 and 1 (2) 
electrons in the up and down spin channel, respectively, for Tb (Ho)). An alternative solution 
to remove spurious effects related to the wrong energy position of Tb 4/ electrons is to treat 
them as valence states using an LSDA+U approach.— Tests performed on the undistorted 
perovskite-like cubic phase using LSDA+U^ on the 4/ states showed that the resulting 
electronic structure is unchanged with respect to the "4/-in-core" treatment. 

The unit cell in the GdFeO"3-like orthorhombic phase shows the Pnma symmetry (20 
atoms/unit-cell, choosing b as the longest axis), with enormous distortions with respect to 
the ideal cubic perovskite (see Fig^): due to the Jahn- Teller instability shown by the Mn 3+ 
ion - with electronic configuration d 4 (t\ g e 1 ) -, oxygen octahedrons are highly distorted and 
tilted (the average Mn-O-Mn angle is ~ 145°, ~ 142° and ~ 144° for TbMn0 3 , HoMn0 3 and 
YMnO"3 , respectively, to be compared with the much larger value of ~ 155° in LaMnO"3 ). 
We recall that for YMnO"3 (and for other R smaller than Tb), the orthorhombic perovskite 
structure is no longer stable, the hexagonal non-perovskite structure competing in stability. 
However, the transition to the metastable orthorhombic Pnma can be obtained by high 
pressure synthesis 15 . The experimental lattice constants and internal positions according 
to neutron diffraction data - very similar for the three compounds, HoMnO"3 YMnO"3 and 
TbMnO"3 - have been used throughout the worki^. In order to separate chemical and struc- 
tural effects, we have also considered: a) LaMn03 , having less distorted octahedrons and b) 
two "artificial" systems, namely i) YMnC>3 in the LaMnC>3 structure and ii) LaMnC>3 in the 
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YM11O3 structure (denoted in the following as [YMnC>3 }i a and [LaMnC>3 }y, respectively). 

III. RESULTS AND DISCUSSION 

Let us first focus on YMnO"3 , TbMnC>3 , HoMnO"3 and LaMnOs in their equilibrium 
structure. In Table U we report the total energies (calculated per formula-unit containing 
one Mn atom) with respect to the FM phase. If we restrict to the FM and AFM - A, 
C and G orderings - (i.e. we neglect possible next-nearest-neighbors AFM coupling), the 
situation for HoMnO"3 , TbMnC>3 and YMnC>3 is common to the well studied LaMnO^, 
where the most stable magnetic configuration, as experimentally observed^ and theoretically 
confirme d 17 ^ 18 , is the AFM-A type 19 . However, the inclusion of the AFM-E in the subset 
of considered magnetic configurations largely changes the scenario: interestingly, among 
the collinear magnetic orderings, i) HoMnO"3 show AFM-E as magnetic ground state, in 
excellent agreement with experiments, ii) both YMnO"3 and TbMnO^ show the AFM-E 
type as stable configuration and Hi) in LaMnO^ the AFM-E is also not too high in energy 
compared to the AFM-A stable configuration. Indeed, this confirms what was proposed 
via model-hamiltonians, i.e. the occurrence of the AFM-E phase adjacent to the A- AFM 
in parameter space and competing with the FM metallic phase as well 8 -. The differences 
between HoMnO"3 , TbMnO^ and YMnC>3 - having very similar structural parameters - 
suggests that, to some extent, the 4/ moment affects the Mn ordering. In order to further 
highlight the effect of the spin-polarization due to the 4/ Tb electrons, we considered a 
configuration in which these electrons are kept in the core, but they are equally shared 
between the majority and minority spins (4 in the up and 4 in the down channels, denoted 
as (4j,4J,)), instead of the "normal" configuration - considered so far - in which there are 7 
and 1 electrons in the up and down spin channels, respectively, denoted as (7|,lj)). The 
results show that the energy ordering of the different phases in TbMnO"3 W' 4 ^ is very similar 
to YMnC>3 , therefore suggesting that the Mn-Mn coupling is appreciably affected by Tb 
4/ electrons. In particular, due to the geometric location of Tb atoms (i.e. in between 
ac-planes), the main differences induced by the treatment of Tb 4/ states occur in the out- 
of-plane interactions, appreciably affected by the 4/ "spin-polarized" or "paramagnetic" 
electronic cloud in the (7 j, 1 [) and (4 j,4 J.) case, respectively. 

At this point, one might infer that the stability of the AFM-E is closely related to 
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distortions (larger in HoMnC>3 , YMJ1O3 and TbMnOs compared to LaMnOa). Indeed, this 
is confirmed by our "artificial" structures: in [YMnC>3 }i a the reduced distortions - compared 
to YM11O3 in equilibrium - restore the stability of the AFM-A (similar to LaMnC>3 ), whereas 
the increased distortions in [LaMn03 ]y - compared to LaMnOs at equilibrium - make 
the AFM-E the ground state. Therefore, the stability of the AFM-E seems to be mainly 
connected with structural (i.e. octahedral distortions) more than with chemical (i.e. identity 
of the rare-earth element) effects. 

In order to investigate the role of electronic correlation on the stability of the different 
magnetic states, we performed some calculations for YMnC>3 in the FM, AFM-A and AFM- 
E, according to an LSDA+U (or GGA+U) approach in the "atomic-limit"— and varying 
the U value from to 8 eV, keeping the U/J ratio fixed and choosing J=0.15*U. In Fig|2] 
a) we show the difference between total energies in the AFM-A and FM as well as in the 
AFM-E and FM. For a complete description of the electronic and magnetic structure, we 
also report the variation of the magnetic moment and of the band-gap as a function of the 
Coulomb parameter (see Fig. E]b) and c), respectively). The bandgap and the magnetic 
moment both increase as a function of U, with a strong and rather weak dependence on 
the Coulomb parameter, respectively. This is expected, since LSDA+U tends to enhance 
the "localization" of Mn d states (thereby increasing the magnetic moment) and to push 
unoccupied states up in energy (thereby increasing the band-gap). As far as the magnetic 
ground state is concerned, our results show that the stability of the AFM-E and AFM-A 
is strongly affected by correlation effects: for values of U larger than 4 (7.5) eV, the FM is 
stabilized with respect to the AFM-E (AFM-A). We remark that the LSDA+U formalism 
might overestimate the tendency towards ferromagnetism, as already noted by Terakura and 
coworkers for LaMnO^ and consistently with our results in YMn03 . In addition, to our 
knowledge and at present, there are no spectroscopic works for YMn03 in the orthorhombic 
phase, on the basis of which one could extract the value of the Coulomb parameter. So, the 
LSDA+U results are mainly shown to warn the reader that the magnetic ground state may 
be affected by correlation effects. We hope our work will stimulate further experimental 
works on distorted manganites to gain additional insights on the accuracy of bare GGA vs 
GGA+U in treating these systems. 

More information on the magnetic interactions between Mn atoms can be gained by 
considering a Heisenberg Hamiltonian 20 and using this model to fit our GGA-calculated 



5 



total energies. In particular, we consider the first- nearest-neighbor (J™) and second nearest- 
neighbor coupling ( J|j inn ) in the ac plane, as well as the first- and second-nearest-neighbour 
coupling out-of-plane (J]_ and Jj_). Before discussing quantitatively our results, let us 
first recall what happens in the well-understood LaMnOs , showing a smaller GdFe03-type 
distortion: the nn FM superexchange (SE) between e g spins competes with the AFM SE 
between ti g spins.— The delicate balance of these two mechanisms crucially depends on the 
MnOg tilting and Jahn- Teller distortions.— Indeed, this is evidenced by our results: for 
the heavily distorted compounds, the FM nn exchange parameters (J|" n ~-l-2 meV) are 
smaller in module^ than the much larger value in LaMnOs (Jjj 1 " ~-9 meV). This can be 
explained considering that, upon large tilting (i.e. in HoMnC>3 , TbMnC>3 and YMnC>3 ), 
the overlap between Mn 3d and O 2p orbitals is not large enough to provide a strongly 
FM SE between Mn cations in the ac plane. In turn, this alters the delicate balance with 
the nnn SE interactions through Mn-O-O-Mn paths along the a direction. Upon significant 
GdFeOs-distortions, the large Ju" nn (of the order of few meV in all compounds) AFM coupling 
prevails, therefore stabilizing the AFM-E ordering. Furthermore, the out-of-plane coupling 
between ti g spins (J]_ of the order of few meV) keeps a strong AFM character along the R 
series, whereas we obtain very small (< 0.5 meV) J\ coupling perpendicular out-of-plane. 
Finally, the octahedral distortions also affect this interplanar coupling (J]_ is smaller in 
LaMnO"3 than in distorted manganites), therefore weakening the assumption, very common 
in model hamiltonians, of restricting the study to a 2D square lattice^, neglecting out-of- 
plane effects. 

Let us now discuss the electronic and magnetic properties of the AFM-E phase, in terms 
of magnetic moments, band structures and density of states. The GGA Mn local magnetic 
moment (~ 3.3 ub) is reduced with respect to the Hund's rule value (4/ig), reflecting the 
hybridization with O atoms (slightly ferromagnetically polarized). The GGA-calculated 
total density of states and band-structure for the AFM-E YMn03 are reported in Fig. 01 
Common to the well studied LaMnOs and shown by the density of states projected on the 
muffin-tin spheres (not reported), there is a strong hybridization between Mn d and O p 
states. In particular, the bands lying just below E F (see highlighted states in Fig|3] (a) and 
(b)) and separated by the rest of the valence band are Mn-0 bands (occupied hybridized 
e 9 -like states). In agreement with the experimentally observed insulating character—, the 
resulting band structure shows an indirect band gap, with the conduction band minimum at 
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r and the valence band maximum at the Brillouin zone edge along the [010] line and - more 
or less degenerately - along the [110]-r line. It is interesting to note that, even without the 
introduction of the on-site Coulomb U repulsion, the AFM-E configuration shows an energy 
gap, in agreement with the experimental observation of an insulating behaviour. When 
comparing with the metallic character of the cubic perovskite-like systems (not shown), 
despite their FM or AFM spin alignments, the effect of the Jahn- Teller distortion and of the 
octahedral tilting is seen to result in a removal of states at the Fermi level by opening band 
gaps.— As shown in FigOc), the precise value of the band-gap is strongly affected by the 
inclusion of correlation effects: for the extreme case of U = 8 eV, the gap increases up to ~ 
2.5 eV. Finally, in order to have some hints on the orbital physics in AFM-E, we plot in Fig. 
El (c) the YMn0 3 total charge density of the occupied states (in the [-0.7;0] eV energy 
range with respect to Ep), which remarkably shows the staggered ordering of the relevant 
e g orbital occurring in the Mn basal plane. 

IV. CONCLUSIONS 

In summary, we have shown via density-functional-based methods that in manganites 
with highly distorted MnOg octahedra, such as HoMnOs , TbMnC>3 and YMn03 , the 
insulating AFM-E phase is stabilized. This is in excellent agreement with experiments, 
where the AFM-E is observed for HoMnC>3 . The calculated exchange parameters show 
that GdFe03-like distortions - that increase along the rare-earth series - reduce the in-plane 
nearest-neighbor FM superexchange interaction along the Mn-O-Mn path. This interaction, 
along with the delicate interplay with the AFM next-nearest-neighbors superexchange along 
planar Mn-O-O-Mn paths is ultimately responsible for the stabilization of the peculiar AFM- 
E spin arrangement. Although the approximations made in this study (such as partially 
inaccurate treatment of correlation effects and atomic relaxations, etc.) should well be 
tested in the near future, we hope that these results will stimulate more works - both from 
theory and experiments - to ascertain the role played by the AFM-E magnetic configuration 
in the manganites framework. 
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Fig. 1 - Picozzi et al. 



FIG. 1: (Color online) (a) The Pnma orthorhombic cell of RMnOs . Black, white and grey spheres 
represent O, Mn and Tb (or Y or La or Ho) atoms, respectively, (b) ac in-plane view of Mn02 
planes, showing the large octahedral distortions and "up-up-down-down" zig-zag spin ordering. 
Solid lines mark the in-plane projected unit-cell. 

spin-spin coupling, which, however, goes beyond the scope of the present work. 
21 I. Solovyev, N. Hamada and K. Terakura, Phys. Rev. B 53, 7158 (1996). 
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TABLE I: Calculated total energy difference (in meV/Mn) with respect to the FM phase. Recall 
that i) type-E consists of intraplane up-up-down-down and interplane AFM coupling; n) type- 
E* consists of intraplane up- up- down- down and interplane FM coupling; Hi) type-A consists of 
interplane AFM and intraplane FM coupling; iv) type-C consists of intraplane AFM and interplane 
FM coupling; v) type-G consists of both interplane and intraplane AFM coupling. Numbers in bold 
denote the ground state. 
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Fig. 2 - Picozzi et al. 

FIG. 2: (Color online) (a) Difference (in meV/formula-unit) between the FM and AFM-E (black 
solid line) and the FM and AFM-A (black dot-dashed line) as a function of the Coulomb parameter 
for YMnG"3 . J is set to 0.15*U. The changes of (b) Mn magnetic moment and (c) band-gap as a 
function of U are also shown for the different phases: FM (dashed line), AFM-A (dot-dashed line) 
and AFM-E (solid line). 
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FIG. 3: (Color online) (a) Band-structure and (b) total DOS for YM11O3 in the AFM-E magnetic 
configurations, (c) Total charge density contour plots for occupied e g states. Arrows denote the 
Mn spins; the solid lines mark the unit cell. 
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